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ABSTRACT 

We consider structured environments modeled by bosonic quantum networks and investigate the probing of their spectral 
density, structure, and topology. We demonstrate how to engineer a desired spectral density by changing the network structure. 
Our results show that the spectral density can be very accurately detected via a locally immersed quantum probe for virtually 
any network configuration. Moreover, we show how the entire network structure can be reconstructed by using a single 
quantum probe. We illustrate our findings presenting examples of spectral densities and topology probing for networks of 
genuine complexity. 


Introduction 

Recently, the study of dissipative or decohering environments as engineered reservoirs has received considerable attention. 
Reservoir engineering aims at turning the environment from an enemy to an ally by modifying the properties of the unavoidable 
noise in order to optimally preserve quantum features such as coherence or entanglement. It is nowadays generally recognised 
that the initial belief of the environment as the major enemy of all quantum technologies must be reconsidered. Theoretical and 
experimental results show indeed that the interaction with the environment, when suitably tailored, can be used to generate or 
protect quantum resources, and holds great potential as a method to increase the feasibility and scale of quantum applications . 1-5 

On the other hand, the extraction of relevant information about a complex quantum system by means of a localised 
quantum probe has also enjoyed growing attention. Besides being generally a less invasive technique, this approach might 
significantly reduce the resources and effort needed to detect, e.g., the complex system’s temperature , 6 communities , 7 state 8 or 
dimensionality . 9 

In this work we apply concepts and techniques of open quantum systems theory to investigate complex bosonic quantum 
networks, thus bridging the gap between these two fields. More precisely, we consider a bosonic probe system locally immersed 
in the network and investigate the control and probing of the network spectral density /(ft)), a quantity embedding both the 
environment structure and the interaction between system and environment . 10, 11 Notably, we demonstrate that not only /(ft)) 
but the entire structure and topology of the network can be obtained by quantum probing, and we give a detailed method to 
accomplish this task. We consider the probe as an open quantum system whose dynamics arises from the interaction with a 
nontrivial quantum environment. 

The study of the spectral density of bosonic quantum systems has been addressed before only for the simple cases of either 
non-interacting oscillators or oscillator chains . 12-16 This leaves the determination and especially control of spectral densities 
of complex quantum networks a largely unexplored area. Previous investigations include using the Landau-Zener transition 
probabilities of a qubit to detect global features of /(ft )), 17 the center of mass motion of a michromechanical oscillator in the 
high temperature limit , 18 the Stokes shift response function of a spin system immersed in a proteomic scaffold and in contact to 
a solvent 19 and the statistics of photons scattering from an open quantum system coupled to a zero temperature reservoir . 20 In 
comparison, our probing scheme has the advantage of being local and of making no assumptions about the temperature or the 
structure of the environment. Networks of spins as well as bosonic/fermionic networks have been also investigated in Refs . 21-23 
in the context of structure probing. However, these methods require a known topology and full state tomography to detect the 
coupling and local field strengths , 21,22 or can probabilistically reveal an unknown topology but are limited to small networks 
only . 2,1 On the contrary, our approach is not based on state tomography, requires minimal post-processing of measured data, is 
in principle deterministic, and can be used for complex quantum networks having hundreds of nodes. 
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Figure 1 . (Color online) Schematics of the network configurations. The white circle is the probe, coupled to a single node 
with a coupling strength k. Couplings between network oscillators are either strong (solid line) or weak (dashed line). Example 
(A) is a periodical chain where every third coupling strength is weaker. Example (B) is a homogeneous chain with one 
short-cut, with all coupling strengths equal. Example (C) is a small-world network constructed from a chain by adding several 
weak shortcuts and example (D) is an Erdos-Renyi random network. 


After introducing the model in detail, we calculate spectral densities of quantum networks of genuine complexity, e.g., 
small-world 24 and Erdos-Renyi 25 configurations, and show how to control some features of their structure, such as the number 
of gaps. Numerical simulations based on exact diagonalization show that the scheme is very accurate for virtually any network 
configuration and temperature. While probing is done in the weak coupling regime, the detected /(©) can be scaled to any 
coupling strength once its shape is known. Finally, we prove that the whole structure and topology of the network can be 
obtained with a single quantum probe by scanning over the network nodes with single and pair-node couplings. 

Results 

The microscopic model. The probe is a quantum harmonic oscillator of unit mass coupled linearly with coupling strength 
k to a node in the quantum network. The probe Hamiltonian, in units of h, is //«, = (/)f + ( 0 ^q 2 s ) /2, where ps and qs are the 
probe momentum and position operators and a>s is the probe frequency. The network plays the role of a finite environment in a 
thermal state of temperature T and consists of N quantum harmonic oscillators of unit mass coupled via springlike couplings, 
each having the same bare frequency ©o- We assume that the probe and the network are initially uncorrelated. 

The Hamiltonian of a generic oscillator network can be given as He = p y p/2 + q r Aq, where p = {p\ .ps- ■■■■Pn) 1 
and q = {<71 , < 72 ,..., < 7 yv} r are the vectors of momentum and position operators and the N x N matrix A has elements A (/ - = 
SjjCOf /2— (1 — Sjj)hij/2, where ft), is the effective frequency of oscillator i, and is the strength of the springlike coupling 
between oscillators i and j. Matrix A is the adjacency matrix (including diagonal terms) of the network and as such completely 
characterizes it. For any network configuration, the network Hamiltonian is quadratic in position and momentum and can 
therefore be diagonalized with an orthogonal transformation K as K y AK = D, where the diagonal matrix D has elements 
Z)„ = £ 2 ?/ 2 , where £ 2 , are the eigenfrequencies of the network. 

By defining new variables P = K y p and Q = K y q, the diagonalization allows us to move into an equivalent eigenmode 
picture of independent oscillators, each interacting with the probe with a coupling strength directly proportional to k. The 
network Hamiltonian in the eigenmode picture is of the form He = Yli=\{P? + £2 2 <2?)/ 2, where /) and Q, are the position 
and momentum operators of the network eigenmodes which have frequencies £2,. The Hamiltonian of the total system is 
H = Hs+He + kll/, where /// is the interaction Hamiltonian. In the eigenmode picture it is of the form /// = qsYj% 1 giQi 
where the dimensionless constants g; describe the coupling strengths between the probe and the eigenmodes and are just 
elements of a row of K. This can be seen by writing the position operator of the node directly coupled to the probe in terms of 
Qi and identifying the weights given by the elements of K with g,. Further details about oscillator network diagonalization can 
be found in, e.g.. Appendix A of reference' 3 by Vasile et al. 

The reduced dynamics can be described exactly by a generalized quantum Langevin equation. In particular, the dissipation 
and memory effects are accounted for by the damping kernel y(t) = ^^ 1 (^ 2 g?/£2 2 )cos(£2,f). The spectral density of the 
environmental couplings is defined as 

/ ( ®) = fl^5 ( ©-£2 ; ). ( 1 ) 

The spectral density and the damping kernel are related as 

rUncix 

J(co) = co y(f)cos(fi)f)dr. (2) 

Jo 
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Figure 2. (color online) The spectral densities (solid line) of the configurations presented in Fig. 1 for a system-network 
dimensionless coupling strength k = 0.01. All networks have N = 200 nodes with a frequency coq = 0.25. The spectral 
densities are probed (circles) with a constant interaction time t = 500. Each circle corresponds to a different value of the 
system frequency, which is varied to sample the underlying J(co). The system is initially in vacuum state and the network 
temperature is T = 5. In the periodic chain (A), strong and weak couplings between network oscillators are 0.1 and 0.06, 
respectively. In the chain with a shortcut (B), all network couplings are set to 0.1. In the small-world network (C), the strong 
and weak couplings are 0.1 and 0.003, respectively. The random network (D) has equal couplings set to 0.05. 


For finite networks, the interaction time between the probe and the network controls the crossover between smooth and 
discrete spectral densities. For t max in eq. (2) smaller than a characteristic recurrence time T/, proportional to the size of the 
network, the spectral density will be a smooth function of frequency. In this regime, finite size effects are avoided and reduced 
dynamics correspond to a probe interacting with a continuum of frequencies, similar to that calculated from eq. (1) in the limit 
N —► °°. In tmax > 1/ regime, finite size effects start to play a prominent role. 

In what follows, eq. (2) will be used to calculate J(co) throughout. We will work in the continuum regime when probing 
and engineering the spectral density, but it will be seen that when probing the topology, one must work in the discrete regime 
instead. 

Engineering of the spectral density. We consider several different network configurations and calculate the respective 
spectral densities. By going beyond simple chain networks a very rich variety of different spectral densities arises. Schematics 
of network configurations considered in this work are presented in Fig. 1, and the corresponding spectral densities are presented 
in Fig. 2. 

A simple homogeneous chain with springlike couplings has a single band. The lowest frequency where J(co) is non¬ 
vanishing coincides with too and the highest frequency increases as the stiffness of the chain is increased, while the magnitude 
is directly proportional to k 2 . Band gaps can be created by periodically tuning some of the couplings in the chain to be weaker, 
which splits the single band to two or more bands. 11 We have checked that the number of bands coincides with the number 
of network oscillators in a strongly coupled group: an example is the trimer chain in Fig. 1 A where several groups of three 
oscillators are weakly coupled as a chain, resulting in three bands as in Fig. 2A. A bigger difference between strong and weak 
network couplings leads to wider bandgaps. 

Instead of tuning many couplings periodically, a similar effect can be achieved by adding just one additional link, or a 
shortcut, to a homogeneous chain, as in Fig. 2B. This leads to a structured J(co) with spikes. The amount of spikes coincides 
with the number of oscillators before the shortcut, while the depth of the pseudogaps between the spikes depends on the 
coupling strength of the shortcut with stronger shortcut leading to deeper pseudogaps. An example with three spikes is shown 
in Fig. 2B. 

We also considered a small-world network and a random network. The former has a single band with a lot of fine structure. 
On the other hand, the random network produces highly structured spectral densities that can change completely from node to 
node. In the example presented in Fig. 2D, the probe is coupled to a node with high connectivity, resulting in a J(co) with a 
wide frequency range. Coupling to a node with low connectivity instead might result in a very narrow J(co). 

The ability to controllably engineer a wide variety of spectral densities highlights that complex quantum networks are 
ideal testbeds for fundamental studies on open quantum systems where, especially in the non-Markovian regime, a number of 
fundamental questions remain still unanswered. 

Probing of the spectral density. We consider the case of a network with an unknown configuration but known temperature 
and show that the local value of the spectral density can be probed by measuring the expectation value of the probe number 
operator (n) = (fl T fl), where a' and a are the system creation and annihilation operators. This is proportional to the average 
energy of the probe. It has been experimentally measured for specific systems such as trapped ions; see, e.g., 27 for a review. In 
the examples, we calculate the evolution of (n) by diagonalizing the network, solving the Heisenberg equations of motion in 
eigenmode picture, and returning to old variables. 

Provided that the network size is not too small, one can easily check by comparison with the exact numerics that, for 
weak probe-network couplings and for times greater than the initial non-Markovian time in the continuum regime t max < T/, 
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Figure 3. (color online) An example of eigenfrequency probing with a homogeneous chain. In (A), the detected spectral 
density is shown for t ss 0.9 t y, where t is the interaction time and t y the recurrence time. The result is a smooth function of 
frequency. In (B), detected spectral density is shown for t ss 3tj. The discrete spectrum emerges and reveals the 
eigenfrequencies marked by dots. In (C), the evolution of ( n{t )) is compared for a>s at an eigenfrequency (squares) and just 2% 
off (circles). The difference is clear and causes the emergence of discrete spectrum in (B). Also depicted is the time Tthermal, 
after which the energy flow is reversed even for a system at an eigenfrequency. The chain, initially in vacuum state, was 50 
oscillators long, with bare frequency ft)o = 0.25 and coupling constant g = 0.1, with a system-network coupling strength 
k = 0.0025. The probe was initially in squeezed vacuum with squeezing parameter r = 1 and phase space angle 0 = n/2. 


the mean value of the probe number operator is well approximated by the expression (n(t)) = e n (/;(())) + N (&)<,-)( 1 — e n ), 
with r = J(a>s)/cOs and where N(a>s) = (e ms/T — l) -1 is the thermal average boson number, T being the temperature of the 
environment, in units of the Boltzman constant k^P The spectral density can be obtained by probing the network scanning 
over different frequencies OJy via the formula 


J{cOs ) = In 



(3) 


where A n{t) = N((0$) — ( n(t )) and t a suitable interaction time. We stress once more that t can be chosen within a wide 
interval, i.e. it should be in the continuum regime and longer than the initial non-Markovian time. A rough estimate of the 
non-Markovian correlation time Tr is given by (O max TR 1 , where (O max is the largest frequency where the spectral density is 
non-vanishing: for weak couplings it is typically very short. 

In Fig. 2 we test the effectiveness of the probing scheme for the spectral densities of the networks presented in Fig. 1 by 
calculating the r.h.s. of eq. (3) with exact numerics for many different values of c o All quantities, namely coupling strengths, 
times and temperatures, are given in arbitrary units as referred to a fixed frequency unit. The probing scheme works well for 
all four examples. In particular it is seen that while the rich structure of Figures 1C and ID would require dense sampling to 
resolve, the samples themselves are accurate. Beyond the examples presented here, we checked the accuracy of the probing for 
a wide range of spectral densities and parameter ranges and found similar results. We also tested probing with different states 
of same energy, finding no differences. 

Probing of the network topology. The spectral density alone does not allow one to determine all properties of the 
underlying quantum network. It can be shown 14-16 that any spectral density that can be produced by a diagonal Hamiltonian 
can be reproduced by a linear harmonic chain with nearest-neighbor couplings only. Since any network configuration, no matter 
how complex, can be presented in the eigenmode picture with a diagonal Hamiltonian, this means that the network shape is 
never unique to a given spectral density. 

In the following we present a general approach allowing for the full characterisation of the network structure and topology. 
Remarkably, our detection method only relies on the use of a single quantum probe and on the measure of the mean excitation 
number (n(t)) at a certain time within a broad characteristic time interval. However, in order to retrieve the complete information 
on the network we need to (i) couple the probe to more than a single node; (ii) measure the excitation number in the discrete 
regime t max > if. The latter point is related to the fact that short-time reduced dynamics are not determined by the entire 
network. 29 This technique is, inevitably, much more expensive than the probing of the spectral density. However, to the best 
of our knowledge, it is the first demonstration that the full information on the network can be extracted by means of a single 
quantum probe. 

We will assume to have some preliminary knowledge on the size of the network, i.e., on the number of nodes N. As 
explained when introducing the microscopic model, a generic oscillator network is determined by matrix A, which provides the 
oscillator frequencies as well as the structure of the network. We will show how to reconstruct both the matrix D and K from 
which matrix A can be obtained via A = KDK 7 . 

To construct D, one needs to detect the eigenfrequencies £2,-, which can be done using an approach similar to J{(£>) probing. 
As explained before, for sufficiently short interaction times, the detected J(co), given by eq. (3), is a smooth function of 
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frequency as if the network had a continuum of eigenfrequencies. With longer interaction times the discrete spectrum emerges. 
This is essentially caused by differences in the heat flow between the probe and the network in the long interaction time regime, 
leading to a significantly larger values of A n(t) when ft )$ is close to an eigenfrequency. The sensitivity of the eigenfrequency 
detection depends on coupling strength k , with weaker values leading to more accurate results. We stress that this behaviour is 
universal to all finite networks. An example is shown in Fig. 3, which also demonstrates the accuracy of the eigenfrequency 
detection when k is weak. Typically eigenfrequency detection is the most robust part of topology probing. 

While £2; do not depend on the node the probe is coupled to, the probe must interact with all eigenmodes to detect all £2,, 
which is not always the case. On the other hand, the couplings to eigenmodes can be changed by coupling the probe to different 
node(s), which allows one in principle to always detect all of them. If £2, are all distinct, as is the case for generic networks, 
knowing N gives the number of £2,, and consequently one knows when all eigenfrequencies have been found. 

The elements of K are detected in two steps. First step is to detect the modulus of K. The coupling strengths to eigenmodes 
are related to the rows of matrix K as gi(j) = Kji, with g, describing the coupling between the probe and the eigenmode i and j 
the index of the node the probe is coupled to. In other words, each node j corresponds to a different Hi and therefore, different 
set of coupling strengths. Thus, the elements Ku can be obtained by probing sequentially each node j of the network at the 
eigenfrequencies £2,, i.e. evaluating decay of (n(t)) at ft)<, = £2,, and using eq. (3) to determine /(£2,). From eq. (1), it can 
be seen that J ( ” ^./ jftjjftxlft) = ^_,k 2 gj. Keeping in mind that we wish to have a value for each of the g t , we approximate the 
integral with a finite sum using J(Q.,) and £2,. By identifying the terms in the summations on both sides and solving for g, we 
arrive at |g,j = -\/2./(£2,')£2;A£2,y (7tk 2 ), where A£2; = £2,- — £2, + i is the sampling interval. 

Second step is to detect the signs of the elements of K by probing pairs of nodes in the same way as previously, e.g. nodes 1 
and 2. In this way one detects Aj,- + Kii\ (with probe frequency £2,) and, by comparing it to the values obtained by the single 
node probing, namely Aj+ \Kii\, one can reveal whether the elements have same or opposite sign. In principle, the choice for 
the tolerance value used here can create some ambiguity, but this can be done after probing is over and the final result remains 
mostly unchanged for a wide range of choices. Running over the eigenfrequencies £2, and running over N — 1 pairs of nodes 
completes the relative signs of the complete matrix. The knowledge of relative signs is sufficient since the adjacency matrix 
A = KDK t is invariant under a global sign change of K. 

Even small errors in the previous steps easily destroy the orthogonality of K, leading to large errors in the reconstruction to 
the overall network structure and topology. To counter this, the orthogonality needs to be enforced. This is straightforward 
to do since it can be shown 26 that any matrix M which has a singular value decomposition M = WNV 1 has a unique nearest 
orthonormal matrix, and it is WV 1 . 

To summarize, one first needs to detect the eigenfrequencies £2,. This determines matrix D. Then, all N nodes are probed at 
eigenfrequencies individually, followed by the probing of N — 1 pairs of nodes, which determines matrix K up to a global sign 
which can be arbitrary. The nearest orthonormal matrix to detected K is calculated and used to calculate the adjacency matrix 
from A = KDK t . The method above can be used for a full reconstruction of a large oscillator network of an unknown topology. 
As an example we show in Fig. 4 a comparison between the original and reconstructed adjacency matrices of a small-world 
network similar to that of Fig. 1C. The accuracy of the reconstruction depends mostly on the magnitude of the adjacency matrix 
elements and the used interaction time t: accuracy improves for larger values of matrix elements and using an interaction time 
in the range %y <t < T thermal, where T t hermai is the time before energy flow is reversed for a system at an eigenfrequency, see 
Fig. 3C. In particular, links that are very weak relative to others can be missed, while a network where all couplings are strong 
leads to best results. The biggest drawback of the method is that the amount of measurements one needs to do scales with N 2 . 
This is a consequence of having an unknown topology with unknown and possibly distinct coupling strengths. Whether it is 
possible to do better without making additional assumptions about the network is still very much an open question. 

Discussion 

We have explored bosonic complex networks from two complementary perspectives. We have shown how, having access and 
control of the network structure, we can engineer a wide variety of ad hoc structured environments. We have demonstrated that 
even small changes in the connectivity of the network, such as the introduction of an additional link, may cause a dramatic 
change in the shape of the spectral density. Highly structured spectral densities are found, e.g., in biological systems and, very 
recently, approaches of open quantum systems theory have been used to understand the loss of coherence in these systems. 34,35 
As biological systems are very difficult to test in the laboratories, it is highly desirable to develop ways to simulate their 
dynamics in a controllable way. Possible experimental implementations of the model considered here include optical modes 
in cavities,^ 2 ions in segmented traps, 10 mechanical resonator arrays, 31 and cluster states created from an ultrafast frequency 
comb. 33 While implementing an arbitrary oscillator network is still beyond the reach of experimentalists, simple networks have 
been implemented with all three platforms and the realization of more complex networks is under active research. 

With a shift in perspective to a complementary view point, we have investigated the possibility of revealing unknown 
properties of a complex quantum network by means of a single quantum probe suitably coupled to the network. We have 
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Figure 4. (color online) The detected matrix A (left), compared to the original one (right) for a smallworld network of 60 
nodes. The initial state for the probe is squeezed vacuum with squeezing parameter r = 1 and phase space angle 0 = tt/ 2, 
while the network is initially in vacuum state. The chain couplings have a coupling strength of 0.2 while the seven extra links 
have a coupling strength of 0.1. Bare frequencies are ft)o = 0.25. 


developed and illustrated a method for efficient detection not only of the spectral density but, in principle, of the entire network 
structure, showing its effectiveness for virtually any network configuration. Our results are the first demonstration of the power 
of quantum probing for complex quantum networks going beyond toy models consisting of tens of nodes. While the present 
results are still a proof-of-principle demonstration, due to the generality of the system considered, they pave the way to the 
study of quantum probing in a variety of physical scenarios. 
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